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We report thermodynamic values of four-point renormalized coupling constant 
calculated by Monte Carlo simulations in the continuum limits of the lattice versions 
of the two-dimensional 0(2) and 0(3) non-linear sigma models. In each case the 
critical index of the coupling constant vanishes, which leads to hyperscaling ( non- 
triviality) . 

A general problem of (Euclidean) quantum field theory (QFT) is whether the n-point 
(n > 2) correlation functions, which are constructed so as to satisfy Osterwalder-Schrader 
axioms Q, are those of the Gaussian model in the continuum limit. If so this QFT is 
non- interacting, and is referred to be a trivial theory 0. It is rigorously proven that one 
cannot construct an interacting continuum limit of the A0 4 theory in the symmetric phase 
in dimensions larger than 4 (D > 4) ||, and is widely conjectured that this is the case even 
in four dimensions Q . A trivial theory, which might be interacting only when a finite cutoff 
is imposed in the theory, cannot be a candidate of a genuine QFT, so it is very important 
to find out whether a theory is trivial or not. Without any proof, it is generally believed 
that a QFT which is not asymptotic free in 4D would be trivial. An important example of 
such a theory is the 4D non-compact QED, and the issue of triviality in this theory remains 
unresolved with variant conclusions of studies. 



As far as the issue of the triviality is concerned, studying the scaling behavior of the four- 
point renormalized coupling constant (g^) is primarily important. A rigorous theorem || 
states that vanishing of is sufficient for triviality in the Ising like ferromagnetic systems, 
and it is conjectured that the same holds for other type of ferromagnetic lattice models [Q. 
While much (numerical) studies on the scaling behavior of g# have been done for the A0 4 
theories in various dimensions || @ and for the Yukawa type interaction ||, little has been 
studied for the 2D O(N) models for N > 2 ||, which share many similarities with 4D lattice 
gauge theories. In particular, the 2D 0(2) non-linear sigma model is not asymptotic free in 
the context of perturbation theory, while 0(iV > 3) non-linear sigma models are asymptotic 
free [|K] . 



In this article, we investigate the scaling behaviors of g^> for the 2D 0(N) (N=2,3) non- 
linear sigma models on a square lattice of linear size L with periodic boundary condition, 
based on Monte Carlo simulations. The model is defined by the action 

A = -/3^a r a J , (1) 

<i,j> 

where (3 is the inverse temperature (/3 = 1/T), and <7j is an N-dimensional unit vector at 
lattice site i. We report compelling numerical evidences that the continuum limits of both 
models are not trivial, irrespective of the property of perturbative asymptotic freedom. 

The four- point renormalized coupling constant (g^) which is basically the connected 
four-point function defined at zero-momentum, and which is constructed so as to be inde- 
pendent of rescaling of field strength, can be written in a system with translation symmetry 



as 



g^ = lim gf = lim U L (UL) D , (2) 
Here, g^ and £^ are respectively the renormalized coupling constant and correlation length 



TT| defined on a finite lattice of linear size L, and Ul is a modified (N component) Binder's 
cumulant ratio defined as U L = [(1+2/iV) < S 2 > 2 -) < S 4 >}/ < S 2 > 2 , with S 2 = | £ 4 a { \ 2 
0- 



For a system displaying power-law type scaling behavior, g^> scales as g^'(t) ~ 
£-2A+7+d^ w ]2 ere the notations are standard (t, for example, is the deviation of the di- 
mensionless temperature from the coupling defined as t = (/3 C — /3)//3 c ). The inequality, 



2 A < 7 + Dv holds |13[ in general for ferromagnetic systems with power-law singulari- 
ties. For the A0 4 models, the inequality has been rigorously proven for D > 4 so that 
lim^o fi^ — > 0, whereas the equality holds for D < 4 (hyperscaling) ||. In 4D, it is rigor- 
ously proven under some mild assumptions 0] that the scaling behavior of g^(t) has such 
a multiplicative logarithmic correction that lim^o <? (£) ~~ > 0. 

For the 2D O(N) (N=2,3) models, we have calculated g^ by employing Wolff's single 
cluster algorithm |L4j] on square lattice with periodic boundary condition. For each model, 
our continuum limit is achieved by adjusting the value of (3 (or, T) so that the correlation 
length in unit of lattice spacing starts to be diverging; at the same time, the lattice spacing 
scales with (3 so that the correlation length multiplied by the lattice spacing remains a 
constant. In general the auto-correlation time for g^ is much larger than that of x or £, so 
much more computational efforts are required for the precise determination of g^' . For a 
given (3 and L, 20-80 different bins were obtained for our calculations, with each bin being 
composed of 10 000 measurements each of which was separated by 8-15 consecutive one 
cluster updating. The statistical errors were estimated by jack-knife method. 

In order to monitor the effect of finite size in the measurements of g^ A \ for each model 
at an arbitrary temperature in the scaling region, we measured g^ by varying L by 10 from 
L= 10 until g{^ did not vary with further increasing of L. For each model, we conclude that 
(as well as \ and £) converges to its thermodynamic value (within the statistical error 
of very precise Monte Carlo data) on the condition that L/£oo > 7 (Figure(l)). According 
to the theory of finite size scaling the following relation holds for the renormalized 
four-point coupling: 

g?\p)=g^(f3)f g (s), s = L/(-M, (3) 
with f g representing a scaling function characterized by g^\ f g (s) has no explicit tempera- 



ture dependence so that the thermodynamic condition holds for any temperature. 

We thus chose L such that L/£oo — 7 for the direct measurement of the thermodynamic 
g {A) (T=1.19, 1.10, 1.04, and 1.02 for the 0(2) model; (3 = 1.5,1.6, and 1.7 for the 0(3) 
model). For the temperatures where the corresponding £ becomes very large, we fix L/£oo 
to be a value much smaller than 7 so that, according to Eq.([3]), thus obtained is exactly 
proportional to its corresponding thermodynamic value. When the corresponding £oo(/3) is 
known at a certain f3, this procedure enables one to obtain accurate thermodynamic values of 
other physical quantities without using sufficiently large L required in direct measurements 



| 15| . Let us illustrate this for the 2D 0(2) model, where very accurate values of £oo are already 
available up to T = 0.98 pf. At T = 1.02 (^ = 26.20(20)), we obtained = 8.86(13) 
using L= 184 and g^ = 4.79(2), so f g (s) is calculated to be 0.541(10) for s = L/^ = 
2.10(2). At T = 0.98 (f ~ 70.5(7)), choosing L=148 makes the value of s the same as 
at T = 1.02 within the statistical errors, so with our Monte Carlo data, g[% = 4.81(2), we 
extract g {4) = 8.89(20) for T = 0.98. For the 0(3) model, we fixed the value of s = 1.165(21) 
by choosing L=40, 75, and 142 for (3 = 1.7, 1.8, and 1.9 respectively. Our final results are 
summerized in Table (1). 

For such broad ranges of correlation length (5.01(2) < £ < 70.5(9) for the 0(2); 11.1(1) < 
£ < 121.9(7) for the 0(3)) we observe that g^\t) ~ t°. These behaviors are qualitatively 
the same as those in the 2D and 3D Ising models [[/J . We thus conclude that the continuum 
limits are non-trivial for both models. In particular, the perturbative property of the non- 
asymptotic freedom in the 2D 0(2) model and that of the asymptotic freedom in the 2D 
0(3) model do not affect the non-trivialities in these models. Especially, because of the 
asymptotic freedom in the 2D 0(3) model, the value of the renormalized coupling constant 
defined at any other value of the momentum must be smaller than our value, g^ = 6.6(1). 
Another remarkable conclusion of the current study is the existence of hyperscaling in these 
models. Hyperscaling relations cannot be defined in our models due to the exponential 
singularities in the critical behaviors of £ and xi nevertheless, the fundamental claim of the 
hyperscaling in the sense that the thermodynamic correlation length is the only relevant 



scale in the scaling region, still holds in two dimensions irrespective of the type of critical 
behavior. 

The author would like to thank Chuck Buchanan for hospitality; Ghi-Ryang Shin and 
Jaeshin Lee for their support in computing. The Monte Carlo simulations were carried out 
on Convex C240 at the University of Arizona. 
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TABLES 

TABLE I. for the 2D 0(2) and 0(3) non-linear sigma models. For the 0(3) model, we 
fixed s = L/^oq = 1.165 by choosing L=40, 75, and 142 for (5 = 1.7, 1.8, and 1.9 respectively. The 
corresponding values of are 2.56(1), 2.56(1), and 2.59(2) respectively, showing the constancy 



of g^' in 


this ranj 


;e of f3 as well. f g (s) 


is calculated to b 


e 0.394(5). 


^((3) =64.6(5) and 121.7(8) 


for I3=l.i 


1 and 1.9 
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0(2) 


T 


1.19 


1.10 


1.04 


1.02 


0.98 
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8.70(16) 


8.80(12) 


8.71(10) 


8.86(13) 


8.89(20) 


0(3) 


(3 


1.5 


1.6 


1.7 
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1.9 




5 (4) 


6.58(13) 


6.70(15) 


6.50(6) 


6.50(11) 


6.57(13) 
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Figure Caption 

Figure(l): / g {A) as a function of at T = 1.10 = 9.32(2)) for the 2D 0(2) model, 
and at (3 = 1.5 = 11.05(2)) for the 2D 0(3) model. The range of L is [10,80] for the 
0(2), and is [10,100] for the 0(3), varied by 10. 
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